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Abstract: We develop an approach for constructing the Baxter Q— operators for generic 
sl{N) spin chains. The key element of our approach is the possibility to represent a solution 
of the the Yang Baxter equation in the factorizcd form. We prove that such a representation 
holds for a generic sl{N) invariant 7^— operator and find the explicit expression for the 
factorizing operators. Taking trace of monodromy matrices constructed of the factorizing 
operators one defines a family of commuting (Baxter) operators on the quantum space of 
the model. We show that a generic transfer matrix factorizes into the product of N Baxter 
Q— operators and discuss an application of this representation for a derivation of functional 
relations for transfer matrices. 
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1 Introduction 



The notion of an 7^— matrix plays a key role in the theory of lattice integrable systems. It is 
defined as a solution of the Yang-Baxter equation (YBE) [1] . The Quantum Inverse Scattering 
Method (QISM) [2-5] allows one to relate an exactly solvable model with each solution of the 
YBE and provides the methods of its analysis. These methods are the Algebraic Bethe Ansatz 
(ABA) [2,6], the method of Baxter Q— operators [7] and Separation of Variables (SoV) [4]. 
The most widespread and well studied of them, ABA, depends crucially on the existence of a 
pseudovacuum state in the Hilbert space of a model. This requirement holds for a majority of 
the models which found applications in statistical mechanics and quantum field theory. Systems 
like the Toda chain [8] or noncompact spin magnets [9], for which the ABA method is not 
applicable, can be analyzed with the method of Baxter Q— operators and SoV. Unfortunately, 
these methods are less developed in comparison with ABA and their application was so far 
restricted to models with low rank symmetry groups. The main obstacle for the effective use 
of the method of Baxter Q— operators is the absence of a regular method of their construction. 
Nevertheless, in recent time there was a certain progress in the field. The Baxter operators 
are known now for a number of models, see Refs. [9-31]. Special attention was paid to the 
analysis of the spin chain models with sl{n) or Uq{sl{n)) symmetry and their supersymmetric 
extensions. In the approach developed by V.Bazhanov, S.Lukyanov and A.Zamolodchikov [11] 
the Baxter Q— operators are constructed as a trace of the special monodromy matrix, the 
auxiliary space being an infinite dimensional representation of the g— oscillator algebra. The 
method was applied to the models associated with the affine quantum algebras Uq{sl{2)) [11], 
Uq{sl{3)) [19] and Uq{sl{2\l)) [31]. The construction relies on the explicit solution for the 
universal 7^— matrix [32-34]. Since the latter is rather complicated, it causes technical problems 
if one attempts to extend the analysis to a general, Uq{sl{N)), case [30]. 

The spin chain models associated with sl{N = 2,3) and sZ(2|l) algebras were analyzed in 
Refs. [28, 35-37] by a different method. It results in much the same functional relations among 
the Baxter Q— operators and transfer matrices as for the (7— deformed models. This method 
relies heavily upon the special representation for the 7^— operator. Namely, it was shown 
in [38] that the sl{N = 2,3) and s/(2|l) invariant TZ matrices on a tensor product of generic 
representations can be represented as a product of two {sl{2)) and three (s/(3) and sZ(2|l)) 
simpler operators. These operators possess a number of remarkable properties which are easily 
translated into the properties of the corresponding transfer matrices. There is no doubt that 
the factorized representation for 7^— operator exists for a general N. The factorizing operators 
were constructed in the explicit form in Ref. [29] for an invariant 7?.— operator acting on the 
tensor product of the principal continuous series representations of SL{N,C) group. The aim 
of the present paper is to show that the factorization holds for an invariant 7?.— operator on 
the tensor product of generic (infinite dimensional) highest weight representations of the sl{N) 
algebra (Verma modules). 

The Verma module over sl{N) algebra can be realized as a vector space of polynomials of 
N{N — l)/2 variables of an arbitrary degree, see sect.O We will use this realization throughout 
the paper. The defining equations for the factorizing operators on Verma modules are too 
complicated to be solved directly for a general N. To find a solution we will use the results of 
Ref. [29] where the invariant TZ operator acting on the tensor product of the principal series 
representation of SL{N, C) group was constructed. This construction uses the properties of the 
intertwining operators for the principal series representations and naturally gives rise to the 



2 



factorized form of the 7^— operator. The building blocks for the 7^— operator are simple integral 
operators defined on functions from L?'{Z x Z), where Z is the group of lower unitriangular 
matrices. We try to interpret these operators as operators on Verma modules. Below we show 
that such interpretation is possible and find explicit expressions for the factorising operators. 

As an illustration let us briefly consider the simplest case of the SL{2, C) group. The 
SL{2,C) invariant 7^— operator acts on the tensor product of two unitary principal series rep- 
resentations, T(^I'^i) (g)T(*2,^2)^ 

[T(^''\g)f]{z) = {cz + dr'%-c-z + d)-'' f , 
s = 1 — s*, 2{s — s) = n, / G L^(C). It can be represented in the factorized form [29] 

Tl{u -V)= V{ui - U2)S{vi - Ul)V{vi - U2)U{UI - V2)S{V2 - U2)U{U2 - Ui) . (1.1) 

All operators U, V, S depend on two spectral parameters, holomorphic and antiholomorphic 
ones, the latter (antiholomorphic) is not displayed explicitly, i.e. U{X) = U{X, A), with A— A = n 
being integer. In Eq. (jl.lj) we have put ui = u + l — si, U2 = u + si, vi = V + I — S2, V2 = v + S2, 
and similarly for the barred parameters, ui = u + 1 — si, etc. The operators U and V are 
integral operators 

[u{x)f]{z,w) =A{x) J d^cr'-'^r'-'^fiz-cw), 

[ViX)f]{z,w) =AiX) j d^r,-^-"^?)-^-'^ f{z,w-rj), (1.2) 

(the factor A{X) is defined in Eq. (j2.36p ) while the operator S{X) is a multiplication operator 

[5(A)/](z, w) = {z- w)-\z - wr^fiz, w) . (1.3) 

For A* = —A all operators are unitary operators with respect to the standard scalar product 
on L^(C X C). It is clear that an interpretation of the operators U, V and S as operators on 
the product of Verma modules causes a lot of problems and hardly possible at all. However, 
for the products of U, V and S operators, R^^^ = VSV and R^"^^ = USU, such interpretation 
exists. Indeed, the action of the operator i?*^^^ on a test function can be represented in the 
formQ 

[R^^'>f]{z,w)=A{u2-V2) I d'^[l-C]"'-'[C]"'f{az-w)+w,w), (1.4) 

where we put for brevity [^]^ = C'^^'^. Let us explore this integral for the case when f{z, w) is a 
holomorphic polynomial in z and w. It is clear that the result is a polynomial again provided 
that the integrals 

4(a,/3)= I d^^eil-Cr'^i^]-^, (1.5) 

where a = 1 + U2 — V2 and /? = V2 — ui, converge for an arbitrary k in some region of the 
parameters a and /3. Let us assume that a — a = n > and (3 — P = m > 0. The integral (jl.5p 

*The derivation of Eq. (|1.4p makes use of the so-called star-triangle relation [39], see e.g. Refs. [9,29] for 
details. 
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converges in vicinity of the singular points ^ = 0, 1, oo if Re a < 1, Re (3 < 1 and Re a + /5 > 1, 
respectively. These conditions can always be satisfied. Thus Eq. (jl.4p defines an operator 
on the product of Verma modules. It could be checked that this operator coincides with the 
factorizing operator obtained in [38]. Next, the integral (jl.Sp is an analytic function of a,f3, 



hence one can extend the domain of definition of operator ()1.4p to arbitrary a, (3. Thus we 
have constructed the factorizing operator on the tensor product of s/(2) Verma modules starting 
from the solution for the SL{2, C) case. 

In what follows we extend these arguments to a general case of sl{N) invariant TZ operator. 
Instead of a direct calculation of integrals (which becomes too cumbersome) we accept another 
approach. An operator on Verma module is determined by its matrix in some basis. However, 
such a description is not very convenient. It is preferable to describe an operator by its kernel 
which is defined as follows, A{z,w) = Ylikn^k{^)-^knh n(^)) where 6^(2) and {w) are basis 
vectors in Verma module and its dual space. The function A(z, w) depends on the functional 
realization of the dual space, but as the latter is fixed there is a one to one correspondence 
between an operator and its kernel. Proceeding in this way we derive the defining equation for 
the kernel of the factorizing operator. Instead of solving this equation we show that the "kernel" 
of SL{N,C) operator calculated in some specific basis satisfies the same equation and, hence, 
provides the kernel of the factorizing operator on Verma module. The factorizing operators in 
the SL(N,C) case possess a number of remarkable properties. The heuristic arguments given 
above imply that the operators obtained by ^^SL{N,C) — > sl{N) reduction" should inherit all 
these properties. 

Taking the trace of monodromy matrices constructed from the factorizing operators one 
defines a family of commuting operators (Baxter Q— operators) on the quantum space of the 
model. The trace is taken over an infinite dimensional space and, in a general, diverges for the 
model with unbroken sl{N) symmetry. The finiteness of the traces can be provided by intro- 
ducing a regulating factor [15, 19, 24, 25] which breaks sl{N) symmetry down to its diagonal 
subgroup. We prove that in this case the traces for the Baxter operators converge absolutely. 
Moreover, the Baxter operators can be identified with the generic sl{N) transfer matrix with 
a specially chosen auxiliary space. Using the properties of the factorizing operators we show 
that a generic transfer matrix can be represented as a product of Baxter Q— operators. This 
representation is quite helpful for the study of the functional relations among the transfer ma- 
trices. Though we will mainly discuss the homogeneous spin chains the approach is applicable 
to the analysis of the inhomogeneous spin chain models. 

The paper is organized as follows: In Section [2] we recall some basic facts about repre- 
sentations of SL{N,C) group and fix the notations. We give also the summary of the results 
of Ref. [29] which will be necessary for further discussion. In Section [3] the sl{N) invariant 
7^— operator on the tensor product of two Verma modules is constructed. To handle some 
technical problems we introduce an invariant bilinear form on a Verma module and describe 
an operator by its kernel with respect to this form. Using this technique we find the explicit 
form of the factorizing operators in the sl{N) case. We prove that these operators obey certain 
commutation relations. In Section H] we construct the Baxter operators and study their prop- 
erties. A generic invariant transfer matrix is defined as the trace of the monodromy matrix 
over an infinite-dimensional auxiliary space. To ensure convergence of the trace we introduce 
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a boundary operator which breaks the sl{N) symmetry of the transfer matrix to its diagonal 
subalgebra. We prove that the corresponding trace over an infinite-dimensional space exists 
and show that a generic transfer matrix factorizes into a product of Baxter Q— operators. 
Section [5] contains concluding remarks. The derivation of some technically involved results is 
presented in the Appendices. 



2 SL{N,C) invariant 7?.-operator: principal series representa- 
tions 

To make further discussion self-contained we give here a brief review of the principal series 
representations of the SL{N, C) group, and then formulate the results of the Ref. [29] which 
will be used in subsequent analysis. 



2.1 Principal series representation of SL{N,C) 

The unitary principal series representations of the group SL{N, C) can be realized on the space 
of functions on the group of lower triangular N x N matrices [40,41]. Namely, let Z_(Z_|_) and 
H^{H-) be the groups of lower (upper) unitriangular matrices and upper (lower) triangular 
matrices with unit determinant, respectively. 



/I 

Z21 1 
Z31 Z32 1 








\Zm ZN2 ■ ■ ■ zn,n-i 1/ 



//ill hi2 hi3 

/l22 h23 
/l33 

V 



hiN \ 

h2N 

hN,N/ 



Almost any matrix g = SL{N,C) admits the Gauss decomposition g = zh. The element 
zi G satisfying the condition ■ z = z\ ■ h will be denoted by zg, so that g~^z = zg ■ h. 
The homomorphism g — > T'^{g), 



[T^{g) mz)=cx{h-^)^z-g). 



(2.1) 



defines a principal series representation of the group SL{N, C) on a suitable space of functions 
on the group Z_, $(z) = $(221, Z2i; -^ai, ^31, • • • , -^A^Af-ii ^A^Af-i) [40,41]. The function a in 
Eq. (12. ip is the character of the group , 



N 



a(/i) = J] /ij 



"-kk 



(2.2) 



fe=i 



Here h^k = {hkk)* is the complex conjugate of h^k-, whereas in general 7^ Uk- We put 
a = (o"!, . . . , (Tat) and will sometimes use notation T*^ instead of T". Since det/i = 1 the 
function OL{h) depends only on the differences a^^k+i = ^k ~ <^k+i and can be rewritten in the 
form 



N-l 



N-l 



cx{h) = \{{^,{h)Y- 



'Ak{h)) 



1 — CTfe^fc+l 



o"fc,fe+l) 



(2.3) 



k=l 



k=l 
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where = a^^k+i — k = 0, . . . , N — 1 are integer numbers, G Z 0. The function 

Afc(M) is defined by 

Afe(M)=detMfc, (2.4) 

where the k x k matrix Mk is fc-th main minor of the matrix M, and Afc(M) = (Afc(M))*. 
That is Afc(/i) = nj=i ^ ^ -^+- We will assume that the parameters 0"^ satisfy the 

restriction 

ai + a2 + ... + aN = N{N- l)/2. (2.5) 
The operator T°'{g) is a unitary operator on the Hilbert space L?'{Z^), 

{^i\^2)= j Dz{^i{z)y^2{z), Dz= J] (fzki, 

l<i<k<N 

if the character a'{h) = a{h) IlfcLi l^fcfeP'^ is a unitary one, i.e. \a'\ = 1. This condition holds 
if ^fc,fc+i + = for A; = 1, . . . , - 1, i.e. 

o-k,k+i = + , o-k,k+i = ^ + , k = 1,2,. . . ,N - 1 , (2.6) 

where nk is integer and is real. The unitary principal series representation is irreducible. 
Two representations and T"' are unitary equivalent if and only if the corresponding param- 
eters (o"i, . . . , ctat) and {a[, . . . , a'j^) are related by a permutation, see for details Refs. [40,41]. 

We will also need the explicit expression for the generators of infinitesimal SL{N, C) trans- 
formations, which are defined in the standard way 

[T"(]l + Y,,, f ] (Z) = Hz) + Y.,k + ^'^) "^(^^ + • 

Here £ik, (1 < ^ < N), are the generators in the fundamental representation of the SL{N, C) 
group, 

(^ifc)nm — ^in^km '^^ik^nm- (^•''') 

The generators Eik{Eik) are linear diff'erential operators in the variables Zmn{zmn)-, (1 < n < 
m < N) which satisfy the commutation relation 

\Eik 1 Enm] — ^kn -^im ^irn ^nk • (^•^) 

The generators Eik admit the following representation [29] 

■^ik — ^ ^ Zkm (^I-^nm ~l~ ^nm '^m^ {z )ni • (^'9) 

Here Dnrm n > m are the generators of the right shifts, ^{z) ^(zzq). 



^From now in, since each variable a comes along with its antiholomorphic twin a we will write down only 
holomorphic variant of equations. 
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Equation (|2.9|) can be written in the matrix form 



E = -z{(T + D)z-^ , (2.11) 

where 

E = ''^Eikeki, , D = ''^Dikeki, cr = ''^akekk (2.12) 

ik i>k k 

and the matrices Cnm {n,m = 1, . . . ,N) form the standard basis in the space Mat(A^ x A^), 
{^nm)ik = ^in^mk ■ The generators D^i satisfy the same commutation relation as -Efci) (Eq. (|2.8p ) 
[Dki , Dnm] = SinDkm — ^kmDni and commute with the generators of left shifts, E'fcj, k > i. 
An explicit expression for the generators of left and right shifts reads 

in Nr. 

\ - a \ - _ a 

Eki = — Z'im jr = Zmk JTZ , (2-13) 

m=l m=k 

' d ^ d 

Dki = - ^ ZimTTZ = ^ Zmk ^ , (2-14) 

m=l m=k 

where Zki = {z^'^)ki and we recall that za = 1. Let us notice here that the operator D^i 
depends on the variables in the A:— th and i—th columns of the matrix z, or on the variables in 
the k—ih. and i— th rows of the inverse matrix z^^ . 



2.2 Coherent states 

In this subsection we describe the system of functions with "good" transformation properties 
with respect to SL{N,C) transformations. This system will play a key role for establishing 
relationship between 7^— operators defined on different spaces of functions. Namely, we define 

N~l 

A'^iz, a) =11 (Afc(az))-i+'^'=''=+i (A,,(az))-^+'^'=''=+i , (2.15) 

k=l 

where Afc(M) = detM^, see Eq. ()2.4p . The function A'^{z, g^^) is nothing else than the 
prefactor a in Eq. (|2.ip . i.e. A*^ is a transformation of unity, 

A-(z,5~^) = T^{g) ■ 1 = a{h-\z,g)) . (2.16) 

We will refer to the function A"'{z, a), with a being an upper triangular matrix, a S Z-^-, as a 
coherent state. For an unitary representation the system of coherent states {A'^(z, a), a G Z^} 
is a complete orthogonal system in L'^(Z). Indeed, it is easy to verify that the integral opera- 
tor D defined as 

[T)ip]{z)= DaA'^{z,a)ip{a), a e Z+, Da = JJ d'^aik (2.17) 

l<i<A,<Ar 

intertwines the unitary representations and T~f, 

T"(5)S) = S)f^(5). (2.18) 
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The representation T'^ is defined on functions on the group 

N-l 



[T'^{g)f]{a)=-f{h)f{a-g), -f{h) = J] (Afc(/i))-^-'^'^^^+i (A,(/i)) . (2.19) 

k=l 



Here for a € Z+ and g E SL{N,C) we put a-g = h-ag, h G It fohows from Eq. (j2.18p that 
the operator DD^ commutes with all operators T°^{g). Since T"^{g) is an operator-irreducible 
representation [40] the operator is proportional to the unit operator. Hence 

j DzA''{z,a) A-(z,a') =Civ(cr) [] 5\aik - a'^^) , (2.20) 
where S'^{z) = 6{x)6{y) for complex z = x -\- iy. For the normalization factor we obtained 



cn{(t) = Yl 



7r2 



l<fc<i<Af ' ' 

As was mentioned above coherent states have good transformation properties. Namely, one 
easily derives 

a-\h+{z, g)) A'^izg, a) = a-\h^{a, g-^)) A'^(z, ag-^) , (2.21) 

where, we recall, g^^ ■ z = zg ■ h^{z, g), a ■ g~^ = h-{a,g~^) ■ ag~^. Equation (|2.2ip can be 
brought into the form 

T^{g) A-(z, a) = f-'^{g~') A-(z, a) , (2.22) 

where the transformation T^"'{g) is given by Eqs. (|2.19|) with the substitution — > — ct/j. 
Thus the coherent state A"'{z,a) satisfies the equation 

=0' (2-23) 



where Eik and Eik, i,k = 1, . . . ,N, are the generators of the sl{N) algebra in the representations 
T"' and T"*^, respectively. 



2.3 Factorized form of TZ operator 

The Yang-Baxter equation (YBE) is an operator equation which is a corner stone of the theory 
of integrable systems. It has the form 

TZuiu - v)TZi^{u - w)TZ23{v - w) = Tl23{v - w)lZi3{u - w)TZi2{u - v) . (2.24) 

The operators act on the tensor product Vi V2 As usual, an operator with the indices ik 
acts nontrivially on the space VjiSiVfc only, (i.e. TZu = Ri2'i^^3 and so on). We are interested in 
SL{N, C) invariant solutions of the YBE. Therefore, it will be assumed that each space carries 
a certain representation of the SL{N, C) group. For a special choice, V3 = V/, where V/ is the 
fundamental representation of SL(N,C), the YBE turns into RLL relation. It has the form 

TZuiu - v)Li{u)L2{v) = L2{v)Li{u)TZi2{u - v) , (2.25) 
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where L{u) is the Lax operator 

L(n) = li + ^ Smn Enm ■ (2.26) 

mn 

With the help of Eq. (j2.9p the Lax operator (|2.26p can be represented as (see Eq. (|2.1ip ) 

L{u) = z{u- a - D)z-^ . (2.27) 

The form of the operators of right shifts does not depend on a representation, hence the Lax 
operator is completely determined by the set of numbers (spectral parameters) {ui, . . . ,un}, 
where = u — cr^., i.e. L{u) = L{ui, . . . ,un). 

A solution of the YBE for the principal series representations of the SL(N, C) group, i.e. an 
TZ operator which acts on the tensor product of the principal series representations H, 
has the factorized form [29] 

TZuiu -v) = PuR^^iui - vi)R^^\u2 - V2) . ..M.[2\uN - vn) . (2.28) 

Here P12 is the permutation operator, Pi2fiz,w) = f{w,z) and the parameters u^, are 
defined as follows, Uk = u — and Vk = v — pk- All operators depend on holomorphic and 
antiholomorphic spectral parameters, i.e. R('')(A) = M('^^(A,A), which are subjected to the 
restriction A — A S Z, but for brevity we display the holomorphic variables only. 
The defining equation for the operator M^^^ is 

M!l^\um - Vm)Li{ui, . . .,Um, ■ ■ ■ Un) L2{vi, . . . , Vm, ■ ■ ■ , Vn) = 

Li(ni, ...,Vm,-- ■,Un) L2{vi, ...,Um,-- . , VAt) M^™^ (li^ - Vm) (2.29) 

and similar for the antiholomorphic Lax operators. Thus the operator M.^^ (um — Vm) exchanges 
the spectral parameters Um and Vm in the Lax operators. It follows from Eq. (|2.29p that M^™) (A) 
intertwines the representations T" and T"™-^ (8) T^™--^, 

M(™)(A)r"(5)®T^(5) = r"™.H5)®2^^'"'-M5)]R^"^(A), (2.30) 

where 

^m,xih) hmm ^mm Ck(^) ; Pm,—xi^) ^mm ^mm Pi^) ■ 

The operator M^"^) is completely determined by the "quantum numbers" of the representations 
they act on, i.e. by the characters cx and /3, (cr and p ) and the spectral parameter A. We 
will display the spectral parameter A as an arg ument of the operator, M(™)(A|q!, /3) M("^)(A), 
and omit the dependence on a and /3, assuming that these parameters are always fixed by a 
representation the operator acts on. 

*We accept the following "standard" parametrization for the characters, a{h) = Y[k=i ^kk''^'' ^kk''~'' ' 
~ Ylk=i^kk''~ ' ^kk*"^ ■ Also, we will usc the letters z and w for arguments of the functions from the 
representation space of T" 8 , f{z,w). 
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The solution of Eq. (|2.29|) can be represented in the form [29] 




(N-l \ /m-1 \ 

Vj(7;j+i - fJUi('Um-Ui) • (2.31) 

j=m y y i=l J 

The operator S(A) is a multiphcation operator 

S(A) fiz, w) = [{w-'z)m]^ fiz, w) , (2.32) 

where [a]^ = a'^a^ . The function [a]"^ is single valued only if A — A G Z. This condition 
is always satisfied in the above construction. The ordered products are defined as follows 
nS^i = AiA2...Am and f^^Ai = A^ . . . Ai . 

The operators Uj(A) = Uj(A, A) and Vj(A) = Vj(A, A) are unitary operators if A* + A = 0. 
They can be expressed in terms of the operators of right shifts, Eq. (j2.13p . acting on z and w 
variables, respectively 

U.(A,A) = (A(:y' (A?m)\ (2-33) 

V.(A,A) = (l^j:l,y (2.34) 
The operators (|2.33|) are well defined and can be represented as integral operators [29], 

[U,(A) <I>](z) =^(A) / d\ [C]-^-^ Hz^) , (2.35) 



where = C C, H = ^ -Ce^+l,ij and 

A{X) = A{X, A) = - r(l + A) /r(- A) . (2.36) 
vr 

The operator U^, with the spectral parameter A = dk^k+i = ffc ~ ffc+ii intertwines the 
representations, and , 

Ufc(afc,fc+i) r"(5) = T-'(5) Uk{ak,k+i) , (2.37) 

where a'{h) = {hkk/hk+i,k+iY''-'"'"^ oi{h). 

2.4 Properties of factorizing operators 

The operators M'-'^^ satisfy a number of remarkable relations [29] 

(2.38a) 
(2.38b) 
(2.38c) 
(2.38d) 

t>M/\ ^ , „ MipH^x ^ , „ \-Tu>i")f\ ^ I ^ ^TtpM^x ^ , „ ^ (2.38e) 





mS?(o) = 


1, 




Mi?(A)Mi?(/x) = 


kS?(a + ^), 




m&)(a)r(?(/.) = 


Mg)(/.)MS^)(A), (n/m). 


m[^^(A)I 


4?(A + m)m(^)(/x) = 


M(f(/.)MS^)(A + /.)M(?(A), 




)^^^^{X-(Jn + Pn) = 


M^a^ (A - CT„ + p„) r!^^ (A - c7„ + Pm 
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These relations are sufficient to prove tliat the 7^— operator (|2.28|) satisfy the YBE [29]. Next, 
taking into account Eqs. (j2.13p . (j2.33p and (j2.3ip one deduces the following commutation 
relations for the operators M^"^^ 

M('")(A) iz-^)kj = {z-^)kj K("^)(A) , for A; > m , (2.39a) 

M(") (A) Wkj = Wkj M^™) (A) , for j <m (2.39b) 

and 

^Si.fcl^^^^W =l^^"'^W^Ki,A:' for k>m+l, (2.40a) 

Dt"+i,k^^""\>^)=^^""\>^)Dt+i,k^ for k<m-l. (2.40b) 

Finally, one can see from the representation ()2.3ip that the operator R(™) depends only on the 
spectral parameters, ui, . . . , Um and Vm, ■ ■ ■ , vn- Namely, it depends on the spectral parameter 
\ = Um-Vm and 

Ui - Um = cTm - (Ti, i < m, and Vm - vj = pj - pm, j > m. (2-41) 

It means that 

M^"^) (A|a, /3) = m(™) (A|a', /3') (2.42) 
if the characters satisfy the following relations 

^ = /(A„(/i), . . . , AM-iih)) and ^ = viA,{h), Am-iih)) . 
That is the ratio a{h)/cx'{h) does not depend on A^, k = 1, . . . ,m — 1 and similarly for /3. 

2.5 Factorizing operators in the coherent states basis 

Let us calculate the action of the operator M^™) in the coherent state basis 

A'^P{z, w\a, 13) = A'^(z, a) AP{w, /3) , (2.43) 
where z,w Z_ and a,/3 € Z^. We will use the following notations 

az = Zad^^aOz , (3 w = w p dyj^p (3^ , (2.44) 

where dz,a-,dw^p are diagonal matrices, Za, wp S Z_ and az-,l3w € -^+- 
Lemma 1. The action of the operator M^™^(A) on the state \2.4^ is given by 

4™)(A)A-''](z,u;|a,/3) =/i?(A) [{Pn. w'' z a~') ^j"" A'^P{z,w\a, P) , (2.45) 
where [a]'^ = a^a?'. The prefactor f^\x) is 



m-l . , \ N 



TT f_nA-A^iK^^ TT ^{Um - Vj) , . . 

n ^ 1) A{uk - u^) ^11^ A{v^ - .,) ' ^'-^^^ 
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where 

A{\) = -i^-^T{l + \)/T{-~X) 

TT 

and 

X = Um-Vm, Um = U-ak, Vk = V - pk ■ (2.47) 

The proof of the Lemma is rather technical and can be found in Appendix lAl Now we want 
to discuss Eq. (j2.45p in more details. First, we notice that the r.h.s. of Eq. ()2.45p is given (up 
to the prefactor f^\x)) by the product of two functions 

JCT^^Jz,w\c^,P) = iP^w'^za^')'^^ (^n(A,(az))'^'=.^+-i(Afc(/?«;))'''=.'=+^-i^ (2.48) 

and 

jCl^ .iz,w\a,P) = [{P^w-'za-%^) ( H (Afc(ztat))<^.,.+i-i(A,(u;t/?t))p-.,.+i-i 



(2.49) 

The function /C (K^) is an (anti)holomorphic function of z, w, a, f3 in the vicinity of the point 
z = 'w = a = P = l. 

Further, it follows from Eq. ()2.29p that the function }C^^ p{z,'w\a, (3) satisfies the following 
equation 

L5'(ni, . . . , . . . Uat) L^{vi, ...,Vm,---,VN) ^™„-,;„,o-,p(^' ^"1"' P) = 

Ll{ui, ...,Vm,...,UN) L2{vi, . . . ,Um, - ■ ■ ,vn) lC^^_^^^„^p{z , w\a, (3) . (2.50) 
Here the Lax operators Li(L2) are given by 

L{u) = li - ^ Cmn Enm , (2-51) 

mn 

where the generators Enm correspond to the representation T~'^{T^p), see Eqs. (I2.22p . ()2.19|) . 
The superscript of the Lax operator (L^, L") indicates the variable it acts on. Equation (j2.50p 
follows directly from Eqs. (j2.23p . (j2.29p . (j2.45p . In the next section we show that the function 
IC^^^ p{z,w\a, P) defines a factorizing operator on the tensor product of Verma modules. 

3 sl{N) invariant TZ operator for generic highest weight repre- 
sentations 

In this section we construct sl{N) invariant solution of the YBE on the tensor product of two 
generic highest weight representations of the sl{N) Lie algebra (Verma modules). 
Let V be a linear space of polynomials of arbitrary degree in z^i, 



V 



^P{z2i,Z3i,...,ZNN-i), deg(P)<oo|. (3.1) 
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The homomorphism 



m<n 

(see Eq. ()2.12p for a definition of Dnm) defines a representation of the sl{N) algebra on the 
space V. The operators E^i are completely determined by the parameters {cJi, . . . , a^}- More 
precisely, they depend on the differences anm = Cn — o'm- In order to stress a similarity with the 
SL{N, C) case, it is convenient to specify the representation by a function a{h) = Y[k=i ^kk~'^'' ' 
where h is an upper triangular matrix with unit determinant. Henceforth, we will use both 
notations, vr" and vr*^, for a representation of the sl{N) algebra. A representation vr*^ is 
irreducible if none of the differences cxjfc = di — ak, i < k, is a positive integer [42,43]. We will 
assume that this condition is fulfilled. 

The highest weight vector, vq, {Ei^VQ = 0, for i > k) oi the Verma module (|3.ip is given by 
vo = 1. For the highest weight A = (Ai, . . . , Aat.i), {{Ekk - -Efc+^fc+i) • vq = XkVo) one finds 
Afc = Tfc+i — cTk + 1. In case if all components of the weight vector A are negative integer, = 
— nfc, > 0, A: = 1, . . . , — 1 the Verma module has a finite dimensional invariant subspace. 
This subspace is a finite dimensional representation of sl{N) algebra which corresponds to the 
Young tableau specified by the partition {ii, . . . where £k = X]i=fc^ ^fc length of 

the k—th row in the Young tableau. 

3.1 Bilinear form and kernel of an operator 

We define the following linear combinations of the Cartan generators E/^^ and the unit operator, 

P N P Q N P g 

Hp = y^(-Efcfc +(7k + k - N) = '^^mk-K = '^^rrik-KZ , (3.3) 

7 1 117 1 Cl^mk ^ 117 1 (yZmk 

k=l m=p+l k=l m,=p+l k=l 

where p = 1, . . . , N — 1 and z = z^^. The space V is a direct sum of the weight subspaces Y^, 
V= Yh = {veY\{Hp-hp)v = 0,p=l,...,N-l}. (3.4) 

Each subspace has a finite dimension. The union of the bases in all Y^ gives a basis in V. 
We will mostly use the following basis 



eniz)=llz^\ (3.5) 



i>k 



where n is a multi-index, n = {n2i, . . . ,nNN-i}- 

Let V be a linear space of polynomials of arbitrary degree in z^i = z'^-, 

Y= \^P{z2i,Z3i,...,ZNN~i), deg(P)<oo} (3.6) 

and if be an antilinear map V ^ V defined by (/9(e„) = = ni>fc ^Ifc*- P^^ = fiYh)- 
Let Q he a bilinear form on the product Y xY such that 

n{v,u) = 0, if veYh and ueYh', h ^ h' (3.7a) 

the matrix fi^m = ^^(en, Cm) is invertible. (3.7b) 
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Let A be a linear operator on the space V and Anm be its matrix in the basis e„, Ae„ = 
em^mn- We wih refer to a function 

A{z,w) = '^en{z){An'^)nmem{w) (3.8) 

nm 

as a kernel of the operator A. The kernel A{z, w) is a (anti)holomorphic function in z{w) in 
the vicinity of the point z = w = 1 {zik = Wik = 0) on condition that the series converges. An 
action of the operator A on an arbitrary vector from V, (which is a polynomial in z) can be 
represented as 

[kP]{z) = Q{A{z,w),P{w)). (3.9) 

The bilinear form is completely determined by a kernel of the unit operator (reproducing 
kernel) 

l{z,w) =^en{z)9.~l^em{w) . (3.10) 

nm 

It is clear that the kernel of an operator does not depend on a choice of the basis. In particularly, 
the kernel A{z, w) of the operator A can be obtained as 

A{z,w) = AI{z,w) . (3.11) 

Let it"' be an irreducible representation of the sl{N) algebra on the vector space V. Hence- 
forth we will assume that the space V is equipped with a bilinear form Qa- such that the 
reproducing kernel has the form 

N-l 

I'^iz, w)=ll (Afc(t«t^))-fc,fe+i-i _ (3_^2) 

k=l 

It is easy to derive from Eqs. (12.11) and (I2.19P that the kernel I"'{z, w) satisfies the equation 

4^)l-(z,t.) = -^lf J-(z,u;), (3.13) 

where Eik = vr'^(eifc) and E^i = 7r~"'{eik) are the holomorphic generators corresponding to the 
representations T'^ (Eq. and T"'^ fEqs. (l2:22]) . (l2l9]l ) It follows from Eq. (f3T3]l that 

the form Qo- (for the irreducible representation) satisfies the properties (|3.7|) . From Eq. (j3.13p 
one derives 

n^{EkiQ, P) + n„{Q, Ek,P) = . (3.14) 

As usual, a bilinear form on the tensor product of two (or more) representations tt"' ir^, is 
defined as 

n„piQl®Q2,Pl(^P2) = n^iQuPl) ■ np{Q2,P2). 

It will be useful to have a more functional definition for the bilinear form Q^-. Let us put 
for any two polynomials P{z) and Q{z) 

B^iQ, P) = C7v(cr) J Dz fi^iz) Q{z) P{z) , (3.15) 



The representation T " was defined on functions on the group Z+, f{a), so that 4f' = Ek^(a)\^= 
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where 



AT-l 

fiUz) = n (Afc(^^^))-'^'=''=+i-^ . (3.16) 

k=l 

Since P and Q are polynomials, the integral in (|3.15p converges in some region of cxfc^fc+i and, as 
will be shown later, defines a meromorphic function of CTj^, which we take for a definition of the 
l.h.s of (jS.lSp for arbitrary aik- Note, that for positive integer ak,k+i the integral (jS.lSp defines 
the invariant SU{N) scalar product. It is easy to check that the bilinear form (j3.15p with the 
measure (|3.16p results in Eq. (j3.14p . Hence the bilinear form B(j, Eq. (j3.15p and the form Qa- 
coincide up to a prefactor. The two forms coincide identically at the following normalization 

c^(^) = 7r-^(^-i)/2 I H a, J (3.17) 

\l<i<k<N J 

(B^iQ, P) = l = n„{Q, P) for Q{z) = P{z) = 1). For example, for = 2 Eq. (I3T5D becomes 
where 2j = ai2 — 1. 

It is useful to extend the space V in order to include into consideration non-polynomial 
functions which are analytic in the vicinity of the point z = 1 [zik = 0). This allows us to 
consider finite transformations of functions from V, P{z) f{z) = a{h) P{zg). The function 
f{z) is no more a polynomial, however, if g is sufficiently close to unity then f{z) is an analytic 
function of z in the vicinity of the point z = 1 { f{z) = Cne„(z) = J2h fh{z), where fh is a 
projection of the vector / to the weight subspace Yh )• For such functions the bilinear form is 
defined as the sum, r2cr(V', /) = J2h ^o-i'tph, fh), in case if the series converges. 

3.2 7^— matrix and factorizing operators. 

Our immediate purpose in this subsection is to construct the factorizing operators IR(™\ m = 
1,. . . ,N which act on the tensor product of two Verma modules V (8) V and solve the RLL 
relation 

M.^j^\um - Vm)Li{ui, . . .,Um, ■ ■ ■ Un) L2{vi, . . . , Vm, ■ ■ ■ , Vn) = 

Li{ui, ...,Vm,-- ■,Un) L2{vi, . . . . . . ,VN)'Ru\'^rn " V.^) . (3.19) 

The parameters Um,Vm are defined by Eq. (j2.47p . It is straightforward to check that the 
operator R2^^''(A) intertwines the representations vr" (E> tt^ (= it"' ® vr'') and vr"™--^ ® vr^™---^ 
(= vr'^' Ovr^'), 

R;™^ (A) vr" ® = tt"-.^ ® M^^) (A) , (3.20) 

where a.^^x{h) = h:^^a{h) and f3„i,~\{h) = h^^P{h). 
"A 



Let TZ^^\z,w\a, (3) be a kernel of the operator m[^^\ 



^i7V](z,u;) = ^7^,p,(7^^^(z,u;|a,/?),V(a,/3)). 
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Using property ()3.14p one easily derives the defining equation for TZ^^\z,w\a, f3) 

Li{ui, ...,Um,... un) L^{vi, . . . ,Vm, - ■ ■ ,vn) 't6^\z, w\a, j3) = 

Ll{ui, ...,Vm,-- ■,un) L^ivi, . . .,Um, ■ ■ . , uat) (z, wj|a, /3) , (3.21) 

wliere A = Um — Vm- We recall here that 

nm 

and similarly for L2 . Let us notice that Eq. (j3.2ip coincides identically with Eq. (j2.50p whose so- 
lution is given by the function /C^^ p{z, w\a^ (3), Eq. (|2.48p . Since the function /C™^ p{z, w\a, /3) 
is an analytic function of z,w,a,P in the vicinity of the point z = w = a = (3 = l 
{zik = . . . = l3ik = 0) it defines an operator on the tensor product of Verma modules V V. 
Thus, we have proven the following statement. 

Lemma 2. The operator M.^^\x): V x V ^ V x V defined by the kernel 

T^^^^piz,w\a,P) = A^ {P^w-^za-')l^I'^{z,a)IP{w,P), (3.22) 
where Am is some constant, solves RLL relation Ii3.19\) . 

As follows from Eq. (j3.22p the operator R^"*^ depends both on the parameters cr and p and 
the spectral parameter A. We will display explicitly the dependence on a spectral parameter 
only. It means that a product of two operators, for instance M.^^\^) ]R.^\x) , written in an 
explicit form, turns into R^™'' {fJ.\cr' , p') M^^^ (A|<t, p), where the parameters cr', p' are determined 
by Eq. I^M)- 

It is easy to see from Eq. (13.22P that matrix elements of the operator R^™^(A) (modulo 
the prefactor Am) are analytic functions of the spectral parameter A. Moreover, they are 
meromorphic functions in cr and p, which specify representations of the sl{N) algebra on 
the space V (8 V. The position of the poles corresponds to the points of reducibility of the 
representations tt"' and tt''. 

The normalization factor Am in ()3.22p is, in general, an arbitrary function of A, a, p: 
Am = Am{X, CT, p). However, the operators Mj^™^ satisfy the relations (|2.38p only if the functions 
Am obey some restrictions. They can be easily read off Eqs. (|2.38p . for instance, 

^„(0,cr,p) = 1, Am{^^,c^' ,p')Am{\,<T,p) = Am{^l + \,<T,p) (3.23) 

and so on. The simplest normalization which satisfy all the requirements is Am = 1. Choosing 
another, SL{N, C) induced normalization, one can get rid of unessential prefactors in some 
formulae for transfer matrices. In this normalization the factor Am reads 

k=l ]=in+l ■' 

where fm{X) = 1 for even m and fm{^) = ^^^^ for odd m. We will assume that the normal- 
ization factor possesses all necessary properties. Its explicit form will be irrelevant for further 
discussion. 
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Let us prove that the operators (X) satisfy the relations (|2.38p . The first of them, 
M.^^\o) = I, follows directly from Eq. (j3.22p . The second one, ()2.38bp . requires a special 
analysis and will be discussed in Appendix [Bj 

Going on to the proof of the relations (|2.38cp . (|2.38dp (|2.38ep we notice that the product 
of the operators M^™) in the l.h.s and r.h.s of the corresponding equations results in the same 
permutations of the spectral parameters Uk, Vk-,Wk in the product of two (three) Lax operators. 
Let us show that two operators which result in the same permutation of the spectral parameters 
in a product of Lax operators coincide up to a normalization. Namely, if 

A{u - v) Li{u) L2{v) =L[{u') L'^iv') A{u - v) , 

B{u - v) Li{u) L2{v) =L[{u') L'^iv') B{u - v) , (3.25) 

then A{X) = cB{\). It follows from Eqs. (f3:25]l that the operator C(A) = B^^{X)A{\) com- 
mutes with the product of Lax operators 

C{u - v) Li{u) L2{v) = Li(n) L2{v) C{u - v) . (3.26) 

We will assume that the operators ^(A) and B{X) and, hence the operator C(A), are analytic 
operators that means that their matrix elements are analytic (meromorphic) functions of A. 
The following lemma states that an operator with such properties is a multiple of the unit 
operator. 

Lemma 3. Let tt^^^ and vr*^^-* be irreducible highest weight representations of the sl{N) algebra 
on the spaces Vi and V2. // an analytic operator A{X) : Vi (8) V2 — > Vi (8) V2 commutes with a 
product of Lax operators 

A{u - v) Li{u) L2{v) = Li{u) L2{v) A{u - v) , (3.27) 

then A{X) = a{X)I. 

Proof It follows from ([STTTj) that 

[^(A),4;)(A)] = KA),4?(A)]=0, (3.28) 

where 

m m 

The space V (8 V can be decomposed into the direct sum of invariant subspaces of the Cartan 
generators, V (8 V = X^/^V/i, dimV/i = < 00. Since 

i^kk +Ekk^Mu)] = 0, the subspace 
Yh is an invariant subspace of the operator A{u). Let e^, n = 1,... ,Nfi be a basis in the 
subspace V/j. Since the representations 7r^^\ vr^^^ are irreducible, the basis vectors are given 
by linear combinations of the vectors 

k>i j>m 

where vq is the highest weight vector, vq = 1. So we write = e^{E^^\E(^^). Now let us 
consider a set of the vectors e^(A) = e'^{E''^\X) , E'-^^X)) . It follows from Eq. (f3?29]) that for a 
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sufficiently large A, the vectors e,^(A) = + 0{1/X), n = 1, . . . ,Nh are linearly independent. 
Hence they form a basis in the subspace Yh- By virtue of Eq. (j3.28p one finds that A{\)e\{X) = 
a(A)e^(A), where a(A) = A{\)vq. Thus, we have proven that 

^(A) = a(A)I (3.30) 

on an arbitrary subspace V/i in some region of A. Due to assumed analyticity Eq. (j3.30p is 
valid for an arbitrary A. Therefore Eq. (|3.30p holds for the whole space Vi (8) V2. □ 



It is clear that the proof of lemma [3 can be easily extended to the case of an arbitrary 
number of Lax operators. Namely, if an analytic operator A(Ai, A2, \m-i) satisfy the equation 

A{\i) Li{u)L2{u + Ai) . . . Lm{u + Xm-i) = Li{u)L2{u + Ai) . . . Lm{u + Xm~i) A{Xi) , (3.31) 

then A{Xi) ~ I. As was explained earlier this result implies that the l.h.s and r.h.s of 
Eqs. (|2.38cp . ()2.38dp and ()2.38ep are equal to each other up to some factor, x. It can be 
easily checked by examining the action of the operators on the highest weight vector, vq = 1, 
that X = 1. We also notice here that Lemma [21 implies the uniqueness of the solution of the 
RLL relation ([339]) . 

The commutation relations ([2.38bp . ()2.40p for the operators M^"*) are vital for our analysis 
of transfer matrices, see Sect. [H The proof of these relations makes use of the explicit form 
of the factorizing operators, Eq. (|3.22p . and the invariance property of the bilinear form fio-. 
First of all, we prove the following lemma: 

Lemma 4. Let g G SL{N) and Ak{g) ^ 0, k = 1, . . . , N - 1. Then 

^a{f-''{g-^)em,en) = n^{em,T''{g)en) ■ (3.32) 

Here en{z) and em{w) = {em{w))* are basis vectors in the spaces Y, V (see Eq. iS. 5]) ). The 
transformations T"'{g),T^'^{g^^) are defined by Eqs. 112.1]) . 112. 19\) . 

Proof. Since Ak{g) 7^ the functions T'^ {g)en{z) and T~'^ {g~^)em{z) are (anti)holomorphic 
functions in the vicinity of the point z = \ [z-ik = 0) 

[T'^iaWMz) = Y,^knek{z), [f"(5-i)e^](z) = B^k ■ (3.33) 

k k 

The matrices A and B satisfy the following relation 

{n-^ B)nrn = {Aa-%m. (3.34) 

which follows immediately from the identity for the reproducing kernel, 

[T^{gm{z,w) = [f^'^ig-')I'^]{z,w), 

(see Eq. (I2.22p ). Inserting (I3.33P into (|3.32p one obtains that the latter reduces to equation 
([33iD . □ 
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Further if both functions Q{z) and [T "'{g '^)Q\{z) are analytic at z = 1 then one obtains 
VL„{f-''{g-^)Q, e„) = Y,<ln ^aicn, T''ig)em) , (3.35) 



where Q{z) = Y.n^nen{z). 

Lemma 5. // the kernel of an operator A has the form 

A{z,a) = r{z,az)1'^ {z,a) , (3.36) 
where the function r{z,a) does not depend on the variables: 

or 

b) {a-^)jn^ j>m 
then 

a) A commutes with z^j, j < m, Az^j = z^j A , 
h) A commutes with zj^ , j > m, A zjf} = zj^ A . 

Proof. The proofs for the cases a) and b) are similar, so we consider the case a) only. Since 
the function r(z, a) depends only on a part of the variables a*^ its expansion in a power series 

r{z,a) ='^rij ei{z) ej{a) (3.37) 

runs only over those basis vectors Cj which lie in the subspaces V/j with the multi-index h = 
(0, . . . , 0, hm, ■ ■ ■ , h]\f-i). Therefore, taking into account Eq. (j3.7p one derives 

n^iriz, a),P{a)) = 17^(r(z, a), P{a)) , (3.38) 

where Hm is a projector to the subspace Ym = Ylh={o o /ijv-i) 

[UmP]ia)=Pia) . (3.39) 

Oikj=0, j<m 

Noticing that the kernel A(z, a) has the form T^"'{z)r(z, a) and making use of Eq. p.35p one 
obtains 



AP{z) = n„i^A{z,a),P{a)j = n^[T-'^{z)r{z,a), P{a) 
Y,njei{z)^^i7^,T''{z-^)P{a)) = ^ r,, e,(z) (i>) , r'^(z-i) P(q)) . (3.40) 

ij ij 

Taking P{z) = ZkjP{z), j < m and noticing that 

r'^(z-i) akjPia) = ZkjUm T'^iz-^) P{a) 
we obtain the necessary result. □ 
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The commutation relations (|2.39p for the operators m[^^(A) is a simple corollary of this 
lemma. The proof of the two remaining relations (j2.40p and ()2.38bp is given in Appendix [Bl 

Finally, with the help of the relations (j2.39p one can check that the operator M(™')(A) 
depends only on a part of the parameters a and p characterizing the Verma modules, namely 



(m) 



(A, (Ti2, . . . , am-lm, Pmm+1, ■ ■ ■ , PN-1n) ■ 



(3.41) 



Proving the properties (|2.38p . (|2.39p . (12.40p we have assumed that the representations of the 
sl{N) algebra the factorizing operators act on are irreducible. This condition can be relaxed. 
Indeed, all relations in question give rise (and are equivalent) to certain equations for matrix 
elements of the factorizing operators M^™^ The matrix elements depend analytically on the 
parameters specifying the representations. Thus these equations hold for arbitrary parameters 
and, hence, for reducible representations. 

On the basis of the obtained results we conclude that the following theorem holds: 

Theorem 1. A sl{N) invariant solution of the YBE on a tensor product of generic highest 
weight representations of the sl{N) algebra can he represented in the factorized form 



45(^2 



V2) . ..r[2\un - vn) ■ 



(3.42) 



Here u-i = u — ai, vi = v — pi. The parameters cr,p specify the representations n^^^ and vr^^^, 
respectively. P12 is a permutation operator. The factorizing operators M.^\x) are given by 
Eq. \3.2S^) and satisfy (at proper normalization) the relations \2. 3^) . 12. 39\) . \2.4-C^ . 



The Verma module is irreducible if none of the differences, cij — fx^, /c > i, is a positive 
integer. In other cases there exists an (in)finite dimensional invariant subspace v, v <ZN. Let 
tt' be a restriction of vr*^ onto the subspace v, n' = tt"'\v, and vr" is a representation induced on 
the factor space W/v. It follows from the RLL relation that the space u V (we assume that 
the representation tt^ is irreducible) is an invariant subspace of the operator 7^i2. Therefore 
the 7^— matrix has a block-triangular form 



^12 (n) 



n[2{u) * 

7^'/2(n) 



(3.43) 



where the diagonal blocks, TZ'i2{u) = 7^i2( 



u] 



and TZ' 



12^ 



iV/v)»Y 



, define new TZ 



matrices on the spaces v®Y and (V/v) ^V, respectively. Thus one can extract 7?.— matrices for 
arbitrary (non-generic) representations of the sl{N) algebra by studying the 7^— matrix ()3.42p 
for reducible generic representations. 



Let us put Um — Vm = A and = v^, k ^ m in Eq. (j3.42p . These constraints mean that 
the characters a and /3 in the tensor product vr" (8) vr^ are related. Namely, one easily finds 
that P{h) = OLjyix{h) = h^'!^a[h). Taking into account that ]r[2^(0) = I one derives 



7eS™^(A) = Pi2Mi7(A) =7^ 



(m) , 



12 



/3=a„ 



(3.44) 



Constructing transfer matrices one refers to spaces in the tensor product V V as a quantum 
(first) and auxiliary (second) spaces, respectively. Clearly, the operator 7?.[^^(A) is completely 
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fixed by the spectral parameter A and the representation 7r"(= vr'^) on the quantum space. 
Making use of Eq. (|3.4ip one gets 

7^(^) (A) = Pi2 (A, ai2 , . . . , a^-i,n , A + a„i , . . . , X + a^-iN) ■ (3.45) 

In the next section we study the properties of the operators defined as the trace of a monodromy 
matrices constructed from the operators 7l^^\x). 

4 Transfer matrices and Baxter Q— operators 

4.1 Transfer matrices 

In this section we discuss properties of transfer matrices for generic sl{N) spin chains. A trans- 
fer matrix is defined as a trace of a monodromy matrix which is given by the product of 
7^— operators. Since we consider infinite-dimensional representations of the sl{N) algebra a 
convergence of the trace is not guaranteed. Except for the N = 2 case, in order to ensure 
the finiteness of traces one has to introduce a regulator (boundary operator) [15, 19, 24, 25] 
which breaks sl{N) symmetry down to its diagonal subgroup. Namely, let us define a new 
7^— operator by 

N-l 

ni2{u,T)=ni2{u,Ti,...,TN-l)= H 4""'^' TZuiu) = T^"' TZuiu) . (4.1) 

p=l 

The operators Hp are defined in (|3.3p and the index 2 { H2 ) refers to the space the operator 
acts on. The operators 7^i2(n, r) obey the YBE 

T^uiu, t) Uisiv, t) n23{v -u)= n23{v - u) 7^l3(^', r) TZuiu, r) . (4.2) 

We define a transfer matrix by 

Tp(n, r) = trp |7^lo(n, r) . . . TIlo{u, t)| . (4.3) 

The index p specifies a representation of sl{N) algebra, 7r''(= vr^), on the auxiliary space. We 
will also label a transfer matrix by a highest weight of the representation in the auxiliary space 
A = (Ai, . . . , Atv-i), Xk = I - Pk + Pk+i, i-e. 

Jx{u, t) = Jp{u, r) = Jp{u, t) . (4.4) 

Let us note that for r < 1 the factor t^^ improves a convergence of the trace since the 
eigenvalues of the operators Hp, p = 1, . . . , N — 1, (see Eq. (j3.3p ) are positive integers. Thus, 
provided that the trace exists, Eq. (j4.3p defines an operator on the tensor product of Verma 
modules Vi <8> . . . (8) V^. The transfer matrices form a commutative family of operators 

[Tp,{u,T),Jp^{v,T)]=0. (4.5) 

It is easy to see that the transfer matrix Tp{u, r) commutes also with the total Cartan gener- 
ators, namely [Tp{u, t), Hp] = 0, p = 1, . . . , N - 1, where H = H^^'^ + + H^^\ 

We will consider homogeneous spin chains only, i.e. assume that the representations of 
sl{N) algebra on the quantum space at each site are equivalent, tt"'^ = tt"'^ = . . . = vr'^^ = ir"'. 

^In Refs. [24] the transfer matrix with an insertion of a boundary operator was defined as foUows T{u, r) = 
tTp{T"°Rlo{u) ...RLoiu)). It is easy to check that these two definitions are essentially the same T{u, r) = 
T^f{u,T^)T-^, where A = J2t=liL - k)Hk. 
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4.2 Baxter Q— operators 

Let us define special transfer matrices, tfie Baxter Q-operators, constructed from the operators 
T^.^™), (fOIj) . Namely, similarly to (fiJ]) we put 

7^ST)(u,r) = r^^7^(^)(n) (4.6) 

and 

Qfc(n + afc,r) =tro{7eS(n,T)...7eg;(n,T)}. (4.7) 

Taking into account Eq. (|3.44p one finds that Qkiu) is given by the transfer matrix for a special 
choice of the auxiliary space 

Qk{u + ak,T)=Jc,^Ju/N,T)=Jx-u(e,-e,^,)iu/N,T), (4.8) 

where the vector is defined by (efc)„ = ^ is the highest weight of the representation 

in the quantum space, = (Jk+i — + 1- 

The Baxter operators Qkiu,T) act on the quantum space of a model and commute with 
each other 

[Qkiu,T),Q^{v,T)]=0. (4.9) 
One easily finds that the Baxter operators satisfy the following normalization condition 

Qfc(fTfc,r) =Pt^, (4.10) 
where H = H^^^ + . . . + H^^^ and V is the operator of cyclic permutation, 

V f{zi,...,ZL) = f{zL,Zi,...,ZL-l). (4.11) 

Let us note also that [Qk{u),V] = [Qkiu),H] = 0. 

Now we are going to prove that Eq. (j4.7p gives rise to a well defined operator on the quantum 
space. To this end it is necessary to show that the trace over an infinite dimensional auxiliary 
space converges. Let Q/fc(n, r) be a monodromy matrix, 

Qfc(n,r) =7^SJ(n,T)...7^^S(n,T), (4.12) 



^ / ^-,m'j...m'^,n' 



Erh,n =en{w) (z''''^) 



k=l 

m'i ...m'r ,n' 



Here we keep the variable w for the auxiliary space and z^'^^ for the quantum space in k—th. 
site. The basis vectors en{z) are defined by Eq. (j3.5p . 

The trace of the monodromy matrix Qfc(ii,r) written explicitely takes the form 

Er^(fc)nm>i,r {fc)imJ,ni . {k)-,m'^nL-i (A-\r>\ 

n=(ni,n2, .■■,"!,) 
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We recall that each summation index, n^, is multi-index, = {(f^fc)^J, i > j}- Let us introduce 
notations, n = (ni, n2, . . . , n/,), and riij = ((ni)ij, (n2)jj, . . . , (ni)jj). First of all, we show 
that all summation indices riij except Uk+ik, ■ ■ ■ ^n^k vary in a finite range, while the indices 

(k) 

Uk+ik, ■ ■ ■ ,i^Nk can be arbitrarily large. To this end we examine an action of the operator TZ^q 
on the basis vector e^; (zi) ^ e^, {w). For brevity we skip the index i.e. zi —>■ mi m, etc. 
Then taking into account Eqs. (j2.39p one finds 






^ n 

i>j>k 


j<k,i 1 





n{z) ® n • (4.14) 

i>j>k J 

From here one concludes that the indices riij for j < k are restricted from above by m^^, 
Uij < m'^p i.e. riij < ''^iij for j < k. 

To prove that rijj is restricted for j > k we use the relation (j2.40p . It takes the form 

= P>f'- (4-15) 

Since for a given m the operators {Dpj^i p)*^*"; where Mp is some number, nullify the vector em{z) 
one derives that {Dp^-^^p)^'^^ Emn{z , w) = for p > k, where Emn(z,w) = Tzl^Q\em{z) (g) e„(ty)). 
The function Emn satisfying these conditions is a polynomial of finite degree (which depends 
on Mp) in Wip, p > k (see Ref. [44], chapter X) hence, riij < M'ij, for j > k. 

Thus we have shown that the summation over Up = {{np)ij}, p = 1, . . . , L in (j4.13p goes in 
a finite range for all {np)ij except {np)k+i,k, ■ ■ ■ , {'np)N,k- Let us also notice that all summation 
indices, {np)ik are of the same order, the difference {np)ik — {nL)ik = {(lp)ik being finite when 
{'nL)ik goes to infinity. It will be shown in Appendix[C]that for the matrix element \JZ^^ {u)]^^' 
in the limit riik — > oo, i = A; + 1, . . . , A^, all other variables being fixed, the following estimate 
holds 

m^l^{u)]t:'\ < c{u)hix\^ ...h%\ (4.16) 

where hj. = '^iLk+i''^ik + 1 and C{u), ak-,----,aN are some constants, which depends on 
a,mij,.... Since t^^ ~ t^'^ ■ ■ ■ t^'^Si the estimate ()4.16p ensures that the series in ()4.13p 
converges absolutely r < 1. Thus we have proven that Eq. ()4.7p provides a definition of an 
operator on the tensor product of Verma modules. 

Let us note that the trace for the operator Qat is given by a finite sum. Therefore, the 
Baxter operator Qn{u-,t) has a finite limit at r ^ 1. The operators QA,.(n, r), k < N, could be 
singular in this limit. However, as it follows from the above discussion, the operator Qk{u,T) 
has a finite limit at Ti ^ 1, i k and Tyt < 1 is fixed. 



4.3 Factorized form of transfer matrix 

The transfer matrix (|4.3p for a chosen quantum space depends on N complex parameters: 
the spectral parameter u and — 1 parameters, — Pk+i, k = 1, . . . , N — 1, specifying 
the representation on auxiliary space. The Baxter Q— operators depend only on a spectral 
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Figure 1: The graphical representation for the matrix element of an operator. 




Figure 2: The graphical representation for the sum Y!>n',m',m" ^t'm''' \.'^^\m' ^nm- 



parameter. The number of independent parameters in the transfer matrix matches the number 
of spectral parameters in a product of N Baxter operators. Below we show that the following 
statement holds 

Theorem 2. The transfer matrix Tp(u,r), t <1, is factorized into the product of the Baxter 
Q— operators 

Tp{u,T) = Qi{u + pi,T) [Vt^Y^ Q2{u + P2,t) ...{Vt^Y^ Qn{u + pn,t), (4.17) 
where H = + ... + 



Proof. The proof of (j4.17p relies on the commutation relations (|2.38cp . Let Tfc(u, r) will be the 
transfer matrix for a special choice of the auxiliary space 

k N-1 

j=i j=k+i 

where cr = (cii, . . . , cttv) are the parameters specifying the representation on the quantum space. 
The transfer matrix Tfc(n, r) can be represented in the form 

Jk{u, t) = tr^^ 7^^;')(n, r) . . . n^l^\u, r) . (4.18) 

The reduced operator TZ'^'^^^ is defined as follows 

n%^^ {u -v,t) = T^°P,oMfo^ {u^-vi)... Rfo^ {uk - Vk) = Pjo {u -v,t), (4.19) 

where uj = u — aj , Vj = v — pj . Let us show that 

Tfc(n,T) =Tfc_i(u,r) (Pr^)"^ Qk{u + Pk,T). (4.20) 

Obviously, the factorization formula (j4.17p is a simple corollary of this result. To prove 



Eq. (|i:20]) let us put A 



u — v) and B 



(k) 

■j = nijo^Ufc - Vk), I.e. 



n%''\u-v,T)=T''°PjoAjBj. 



(4.21) 
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Figure 3: The graphical representation for the reduced transfer matrix T^, the r.h.s of 
Eq. KWf . 
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Figure 4: The graphical representation for the reduced transfer matrix T^, the I.h.s. of 
Eq. dUl]). 

Matrix elements of the operators on both sides in Eq. (j4.20p are given by some sums and one 
has to show that they are equal. To this end it is convenient to use a graphical representation 
for the sums. Let us denote the matrix element A^^' of the operator Aj(Bj) by a box with 
four legs as shown in Fig. [TJ The line connecting two boxes will imply a summation over the 
corresponding index. The operators ( ) will be denoted by an insertion of black (white) 
circle in the corresponding line. As an example, we have given the diagrammatic representation 
for the sum En',m',m" ^"';?;"[r^]™;' B^'^' in Figure [2 

The graphical representations for the r.h.s. and I.h.s. of Eq. (j4.20p are shown in Figs. [3] 
andlU respectively. To obtain the diagram shown in Fig. [3] we have used the commutativity of 
the Baxter operators with the total diagonal generators, H = H^^^ + . . . + H^^^ and represented 
the r.h.s. of Eq. (|i:20]) as Tk-i{u, r) Qk{u + pk, r) t'^ . 

We have already shown that the trace in Eq. (14. 7p converges absolutely for r < 1. Let 
as assume now that the factorization formula ()4.20p holds for the transfer matrices Tk{u,p), 
k = 2, . . . ,p — 1 and that the traces for T^iu, p), k < p — 1 converge absolutely. Then it can be 
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shown that Eq. (|4.2U|) holds for k = p and the corresponding trace converges absolutely. First 
of all, let notice that the summation over indices ki, . . .k^ in 



■ ,mL 



Er ,m^,...,m . 1 -Hik-i,...,kL 

i'fc-lJfcl,...,fei. K QkT \mi..\rr 
k\,...,kL 

goes in a finite range. The blocks Aj and -Bj+i can be interchanged with the help of the 
commutation relation (j2.38cp whose graphical form shown in Fig. [5j Using this identity and 



taking into account that t^^'^^'^'Bjt 




one can transform the sum depicted by the diagram in Fig. [3] into the sum in Fig. [H 

It is clear that the trace for will converge absolutely, if it is the case for T^-i and Qk- 
Since Ti(u, r) = + t) the trace for Ti converges absolutely. Therefore the factorization 
formula (j4.20p and the absolute convergence of the traces for for k > 1 will follow by 
induction over k. This completes the proof of the theorem. □ 



4.4 Fusion relations 

It was long known that the transfer matrices satisfy a set of functional relations which are 
usually referred to as fusion relations. The fusion relations for the compact spin chains were 
thoroughly studied in the literature, see e.g. Refs. [45-51]. Below we show how the represen- 
tation (j4.17p can be used to obtain some functional relations for the transfer matrices. 

Let us consider the product of two transfer matrices, Tp(u, r) and Ttj(z;,r). Both of them 
can be represented in the factorized form (j4.17p . Since all operators in ()4.17p commute with 
each other one can interchanges the Baxter operators Qk{u + pk) and Qk{v + tuk)- The new 
products can be identified as the transfer matrices, e.g. 

(p^W)-(^-l) + -Q + ^ Tp'{u',T) , (4.22) 

where u' = u — 5/N, p'- = pj + 5/N, for j ^ k and p'^ = pk — 5 {1 — l/N), 5 = u — v + pk — uJk- 
Therefore one obtains the following relations 

Tp(n, r) T^ {v,t) = Jp, {u - 5/N, r) {v + d/N, r) , (4.23) 

where tOj = ujj — 5/N, j q and o;^ = Wfc + 5(1 — l/N). Changing the notation Tp(n, r) — > 
Tx{u,t) = Tai,...,Ajv_i (^^1 7") (Tt^(z;,r) — > T^(u,r)) where A (//) is the highest weight in the 
auxiliary space, = Pk+i — /Ofc + 1, one rewrites relation ()4.23p in the form 

'^\{u,t)J^{v,t) = Tx^6^{^^_^-e,){u - h/N,T)J^^s,{e,_,-e,)iv + h/N,T) , (4.24) 
where is — 1 dimensional vector {ek)i = 5ik and 

N-l ^ N-1 

5k = u - V + "^{pp - Xp) + — p{pp - Xp) . (4.25) 

p=k p=l 

The fusion relations ()4.23p . ()4.24p remain valid and for inhomogeneous spin chains. 
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Figure 5: The graphical representation of the permutation identity (|2.38cp . 
4.5 Inhomogeneous spin chains 

Let us explore modifications which appear in a general case of inhomogeneous spin chains with 
impurities (for the sl{2) case see Ref. [35]). The transfer matrix is defined as 

Jp{u, r) = tVpTZioiu + a, r)7^2o(^x + 6, t) . . . 7^Lo(^x + Cl, , (4.26) 

where {^i, • • • are impurity parameters and the quantum space in k — th site carries the 
representation vr^'^) of the sl{N) algebra. 

As in the case of the homogeneous spin chain the transfer matrix (j4.26p can be represented 
in the form (j4.17p . However, properties of the factorizing Q— operators change drastically. First 
of all, let us notice that in the case of an inhomogeneous chain it is not possible to choose an 
auxiliary space representation, ttP, such that the operators TZ^q{u) map -kP (g) tt*^" onto itself, 
simultaneously for all 77. Therefore, the quantum numbers of the auxiliary spaces at different 
sites are different. 

Let us consider the monodromy matrix constructed from the operators TZ^a]p^ {u, t) 

Q,,(n, r) = 7e^%^ (li + Ci, r) . . . u'i^lp^ (u + Cl, r) , (4.27) 

where aj and Pj are the quantum numbers of the quantum and auxiliary spaces at the j— th site 
and Cj are the impurities. In general, the monodromy matrix intertwines the representations 
with different quantum numbers 

Qfc(w, r) vr'' vr'^i (g) . . . (g) vr*^^ = vr'' (g) tt""! (g) . . . vr*^ Qk{u, r), (4.28) 

where we put p = p^. To be precise, for t ^ 1 the above relation holds only for the generators 
from the Cartan subalgebra. Provided that the representations ttP and ttP coincide (p = p), 
the trace of the monodromy matrix Qfc(n, r) = trQ^.(77, r) intertwines the Cartan generators 
of the representations n""^ (g) . . . (g) vr*^^ and tt"'^ (g) . . . (g) vr'^^ . The condition p = p fixes the 
quantum numbers, p„, of all auxiliary spaces in ()4.27p . We recall that 

Rio{u) : tt" ® vr^ ^ 7r"*>" (g vr^*^.— , 

with cxk^u{h) = h^^a{h), jB^-uih) = hlj^P{h) (correspondingly, TZ^iq {u) = Pi2'U.^q{u) maps 
7r"(g)7r^ 7r^'='-"g)7r"'='"). Thus the parameter p„ in Eq. ()4.27p is determined by the characters 
{cin-i)k,u+(n-i- turn the quantum numbers or^, Eq. (j4.28p . correspond to the character 

an{h) = hil-^"+'an+i{h). 
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Let us notice that the character an{h) does not depend on the spectral parameter u and is 
determined by the impurity parameters C only. Since one can always assume that Y2Ck = 
the operator Qfc(n, r) is uniquely determined by two sets of parameters, S = (cri, . . . , <t/,) and 
E = (oTi, . . . ,d-L), 

Qfe(tx,r|C,S) = Qfc(n,T|E,S). (4.29) 

We will not display spins S, S assuming always that the operators are multiplied in a covariant 
way 

Qk{u,T)Qjiv,T) = Qk{u,T\i:,^)Qj{v,T\^,i:) . (4.30) 

Having put the impurities parameters in Eq. (j4.27p to ("j = — {<^j)k we, finally, define the 
operator Q^, as follows 

Qfc(u, r) = t^ll^^lp^ {u + ii- (ai)fc, r) . . . TZ^^lp^ (n + - (fTL)fc, r) , (4.31) 

where the representation ttP^ = ir^" on the auxiliary space at n— th site corresponds to the 
character /3„ = {ctn-i)k,u+^n-i-{crn-i)k' Quite similar to a homogeneous spin chain one can 
show that the transfer matrix for an inhomogeneous spin chain can be represented in the 
factorized form 

Jp{u, r) = r-(^-i)^Qi(7x + put)V-^Q2{u + P2, r)p-^ . . . p-^Qjv(n + pn,t) . (4.32) 

The operators Qfc(u,r) satisfy the following relations d 

Qfc(n, r) Qk{v, r) =Qk{v, r) Qk{u, r) , (4.33) 
Qkiu,r)r-^Qniv,T)=Qn{v,T)V-^Qkiu,r), for k^n, (4.34) 

which follow from the properties of the factorizing operators, Eqs. (j2.38p . These relations allow 
to rearrange the operators in the product of T matrices in an arbitrary order. However, 
the operator alone cannot be considered as a "good" operator on the quantum space of 
the model. Only the product of all operators Qk has the necessary invariance properties. 
Therefore, it is reasonable to identify the Baxter Q^— operator with the transfer matrix (j4.32p . 
where only one operator Q^. depends on a spectral parameter. Namely, let w = {wi, . . . ,W]\f) 
and fi = {fii, . . . , /i7v-i), where pk = Wk+i — Wk + 1- We define the Baxter operator as follows 

=r-(^-i)^Qi(w;i, r)p-i . . . Qk{wk + u, t)V-^ . . . V-^QN{wN,r) . (4.35) 

The definition is not unique in a sense that the Baxter operator depends on the arbitrary 
parameters w (/x). In the case of the homogeneous spin chain there is a distinguished choice, 
w = a, which ensures that Qk{crk) = Vt^ . However, it is clear that other variants are also 
feasible 0. 

"it should be noted that, despite appearance, Eq. (|4.33|l is not a commutation relation. Indeed, if written 

explicitly Eq. (|4:33|l becomes Qfc(M, t|E, E)Qfc(w, rlE, E) = Qfe(i;, r |E, S')Qfe(w, r|E', E). 

** The choice w ^ a could be useful for the analysis of the homogeneous spin chains with a finite dimensional 
quantum space. 
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New operators, (u, r), acting on the quantum space of the model, form a commutative 
operator family. The normalization condition ()4.10p which holds in a homogeneous case is 
replaced now by 

Q'k'\wk,T)=JUO,r)=Z^. (4.36) 

Using the commutation relations (j4.33p it is straightforward to derive that the operators 
(n, r) possess all the properties which hold in a homogeneous case. Namely, the generic 
transfer matrix factorizes into the product of Baxter operators 

Jp{u, r) = Z^-^ qS"^ {u + p,,t)... qP {u + pn,t). (4.37) 

The fusion relation for the generic transfer matrices has exactly the same form as in homoge- 
neous case, see Eq. ()4.24p 



5 Summary 



In this paper we have developed an approach which allows to construct the Baxter Q— operators 
for a generic sl{N) spin chain. We have proven that the sl{N) invariant 7^— operator on a tensor 
product of Verma modules can be represented in factorized form. The factorizing operators 
have an extremely simple form (|3.22p and possess a number of remarkable properties (j2.38p . 
For the homogeneous spin chains we have defined the Baxter Q— operators as the trace of the 
monodromy matrix constructed of the factorizing operators. We have shown that the Baxter 
Q— operators can be identified with the transfer matrices for the special choice of the auxiliary 
space. This definition of the Baxter Q— operators (see Eq. ()4.35p holds for inhomogeneous spin 
chains with impurities as well. 

Many of the properties of the Baxter Q— operators and transfer matrices follow readily 
from the properties of the factorizing operators, Eqs. ()2.38p . (|2.39p . In particularly, we have 
shown that the generic transfer matrix is factorized into the product of N different Baxter 
Q— operators. This representation for the transfer matrix together with the commutativity 
of the Baxter Q— operators results immediately in certain functional relations for transfer 
matrices. 

Another type of fusion relations involves the transfer matrices with a finite dimensional 
auxiliary space. We recall, that Verma module Vp has an invariant finite dimensional submod- 
Pk+i — Pfc + 1 = —nk <Q,k = l,...,N — l. ( This representation corresponds 



ule Vp if Xk 



to the Young tableau specified by the partition {^i, ^2) • • • i ^N-i}^ where £fc = Ylf=k^ '^i is the 
length of the A:— th row in the tableau. ) Let tp(u, r) be a trace of the monodromy matrix over 
such finite dimensional spce 



tp{u, t) = tr„p 7^lo(n, r) . . . 7^Lo(^^, t) 



(5.1) 



Assuming that the normalization of the factorizing operators are chosen according to Eq. (j3.24p 
one can obtain the following determinant representation for the transfer matrix (jS.ip 



tp{u,T) = (Pr^) 



-Af+l 



Qi(u + pi,r) Qi{u + P2,t) ... Qi{u + pn,t) 

Q2{u + Pi,t) Q2{u + P2,t) ... Q2iu + PN,T) 

Qn{u + Pi,t) Qn{u + P2,t) ... QNiu + pn,t) 



(5.2) 
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The proof of (|5.2|) is based on the Berstein-Gel'fand-Gel'fand resolution of the finite-dimensional 
modules and will be given elsewhere [52]. Equation (|5.2p gives rise to a variety of functional 
relations involving the Baxter Q— operators and (in)finite-dimensional transfer matrices. The 
simplest of them are the so-called Wronskian relation and Baxter equation. One easily derives 
from Eq. (j5.2|) the Wronskian relation, which in A notation, tp(u, r) — > tx{u, r) = tAi...Aiv_i i^^ 
reads 

{VT^f^^ to...o{u,T) = det \Qk{u + N -j,T)\k,j=i,...,N: (5.3) 

where the transfer matrix to. ..0(^1 7") is proportional to the unit operator on quantum space. 
Let us put 

to{u,T) = tN{u,T) = to,,,o{u,T), and tk{u,T) = to,,,^i^,„o{u,T) . (5.4) 

That is the transfer matrix tfc(n, r) is given by a trace of a monodromy matrix over a finite di- 
mensional auxiliary space which corresponds to the Young tab leu with one column and fe-rows. 
Following the lines of Ref. [37] one can derive the self-consistency equation (Baxter equation) 
involving the Baxter Q— operators and the finite dimensional transfer matrices, tfc(ti,r). It 
takes the form of the A^— th order difference equation 

N 

Y,i-'^)^tk{u + k/N,T)Qj{u + N-k,T) = 0, (5.5) 

fe=0 

which, due to Eq. (|4.8p . can be considered as a fusion relation involving the finite and infinite 
dimensional transfer matrices of a special type. Equation ()5.5p is a generalization of the 
standard sl{2) T-Q relation. 

Thus the operators Qk{u) possess the following properties 

• Form a commutative family [Qfc('u), Qj(f )] = 

• Commute with all transfer matrices [Qfc('u), Tp(f )] = 

• Satisfy the A^— th order difference equation ()5.5p infolving the finite dimensional transfer 
matrices. 

The operators with such properties are usually referred to as the Baxter Q— operators, which 
justifies using this name in the previous sections. 

Let us notice that relations similar to (j5.2p . (j5.3p . (j5.5p hold for the spin chains with the 
affine Uq{sl{n)) symmetry algebra [19,30,53] (see also [54,55]). 

It follows from Eqs. ()4.17p . (15. 2p that the Baxter operators encode the full information 
about the system. Provided that the eigenvalues of the Baxter operators are known one can 
restore eigenvalues of all transfer matrices. For low-rank symmetry models the knowledge of the 
eigenvalues of the Baxter operators is sufficient to restore the wave function of the system. This 
can be done with the help of the Separation of Variables method developed by E. Sklyanin [4]. 
(Applications for the specific models can be found in Refs. [9,20,23,56-59].) 

The Hamiltonian of the system is defined as a logarithmic derivative of the transfer matrix 
with the same representation in the quantum and auxiliary space. It can be represented as the 
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sum of pair-wise Hamiltonians 



du 

where 



H = — log T^(n, r) = V r^^ Hkk+i r'^'^ , (5.6) 

k=l 



T-Lkk+i = ^fcfc+l(0)^7^fcfc+l(n)|„=o . (5.7) 

Using Eqs. (13.420 and (I3.22p one easily finds the following expression for the kernel of the 
pairwise Hamiltonian (j5.7p 

H{z, w\a 13) = log f n {(iu,w-^za-^)mr^ a)l''{w, (3) . (5.8) 

\m=l / 

Let us notice that for the sl{2) case the argument of the logarithmic function turns into the 
invariant ratio {1 + zP){l + wa)/{{l + za){l + wf])). (Here z = Z21, w = W21 etc.). The 
eigenvalues of the Hamiltonian (|5.6p can be expressed in terms of the eigenvalues of the Baxter 
operators as follows 

L , 

E = y^ — log Qk{u + ak) . (5.9) 

^-^ rill u=0 



du 

k=l 



For a generic representation of the sl{N) algebra it is not possible to define an invariant 
scalar product and the Hamiltonian ()5.6p is not Hermitian in general. Of special interest is a 
situation when the Verma modules possesses an invariant submodule which admits an invariant 
scalar product. It could be a finite dimensional subrepresentation which admits the SU{N) 
invariant scalar product. There are also infinite dimensional invariant subspaces which can be 
equipped with invariant scalar products. They can be identified with unitary representations of 
the noncompact group SU{m, N — m). The spin chain with the SU{2, 2) symmetry group, for 
instance, is relevant for a description of the anomalous dimensions of a certain class of composite 
operators in Quantum Chromodynamics [60]. We hope that the approach developed here will 
be useful for an analysis of this type of models. 
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A Appendix: The SL{N, C) factorizing operator in the coherent 
state basis 

We start the proof of Eq. (j2.45p with a remark that the coherent state (j2.43p can be represented 
in the form 

A''P{z,w\a, (5) = T'^{a-^)Tf^{l3-^) ■ 1 . (A.l) 
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Let us apply the operator M^™) to the function 

<^{z,w) = j DaDl3f{a,(3)A'^P{z,w\a,f3), (A.2) 

where f{a,l3) is a smooth function with a finite support. 
The evaluation of M^™^ <1> is based on the identity 

= 2^«™.A(c,~i)r^m,-A(/3-i)W2 [{w-^(3a-'^z)m]^'Wi, (A.3) 

where 

^jr=m y Y i=l 

i=l / \i=m J 

To derive (|A.3p one uses the following identities, 

WiT'^{a-^)T^{p-^) =r"'(a-i)r^'(/3-i)Wi, (A.6) 
[{w-'z)m]''T"'{a-')T^'{P^') =T^" {a~')T^" {l3~') [{w' pa'' z)m]\ (A.7) 

W2r""(a-^)r^"(/3-^) = r"™.^(a-^)r^'-.-M/3^^)Wi, (a.s) 

where Q:"(/i) = h:^^^a'{h), P"{h) = h^j^l3'{h) and 

m N—1 

cx'ih) = hll- n Kk"-''" ^(h) , (5'{h) = hjf^- n ^^r^' m . (A.9) 

The identities (1A.6|) and (jA.SP follow directly from the intertwining relation (j2.37p . To derive 
()A.7p it is sufficient to notice 

{w~j^^(5a~^Zc,)Ni = {dw,f3 Pw z a^'^ d^}a)Ni = {d^,p)NNiw~^ z)m{d^),)n , (A. 10) 

where we recall that a ■ z = Za dz,a otz- 

Let us now calculate W2 [{w~^ (3 z)n To this end one can use the integral repre- 
sentation for the intertwining operators (|2.35p . Let us show that 

' N-l \ /n-1 \ 

ll^j{X-Pm,j+i)] [{w~'Pa-'z)Ni]^ I n^i^'^-.i+i) I = rm{X)[{w-'Pa-'z)mi]\ 

(A.ll) 



where 



(\\- TT M^^Prnjl_ T-r A{U^ - Vj) 
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and the parameters Uk,Vk, X,ak, Pk are related by (j2.47|) . We represent l.h.s. of (|A.lip in the 
form 

H A{pmj)Aipjm + A) y d\]_^ j 
j=m+l 

where 

w'^ =W (1 - CmCm+lm) ... (1 - ^AT-ieATAr^i) . (A.14) 

It follows from ()A.3P that the integrations in Eq. (jA.lSP have to be carried out in the following 
order • • • , ^A^-i; ^at-I' • • • ' ^'m- Since the integrals converge, we rearrange them and carry 
out the integration over (^at-i, ^^_;^), then over {£,n-2-,S!n-2)-' ^^^^ (Cm,?m)- Each 

integration is reduced to the standard integral 

j d'e I d'i [e' - [e - x,]^ = -j^^^^^ ■ (a.i5) 

(We recall here that [aY" = a^{a*)^.) Then collecting all factors and taking into account that 
{eNN-\ ■ ■ ■ e-m+imW'^ Pa^'^ z)ni = {w~'^l3a~^z)mi one gets Eq. (jA.lip . 

Quite similarly one obtains 

m— 1 \ / m— 1 \ 

[] U,(A - Gi^) \ [(iy-i/3a-iz)„i]^ \ \{ Ui(a,„) j = pU\)[{w-^ [3a-^ z)^^]'' , (A.16) 
where 

m— 1 N m— 1 ./ ^ 

Using these results one obtains from (jA.SP 

=&W [{l^u, z a-')mm]^ A''Piz,w\a, p) . (A.18) 



B Appendix: Proof of the relations ( 12.401) and ( I2.38bl) 



In this Appendix we give the proof of the relations (|2.40p and (I2.38bp for the factorizing 
operators M^™). We start from Eqs. (j2.40p and represent the operator D^+i^k in the form (see 
Eq. (imi)) 

Dk+l,k = —"^^{z ^)ki Zn,k+1 Eni ■ (B.l) 
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Since the relations (|2.4U|) hold for the SL{N,C) factorizing operators one can derive the fol- 
lowing equation for the kernel IC^^^ p{z,w\a, (3), 

in 

where k > m+1 and we have used the commutation relations (|2.39|) for the SL{N, C) operators. 
Taking into account that the kernel of the sl{N) factorizing operator M('")(A) is given by the 
function IC^^^ p{z,w\a, jd) and making use of Eqs. ()3.14p one can easily check that Eq. (|B.2p 
gives rise to the following equation 

41.]K^"^^(A)=K(™)(A)Z?il, (B.3) 

for the sl{N) factorizing operator M*^™'^(A). 

The relation ()2.38bp is equivalent to the following equation for the kernels 

This equation is a consequence of the commutation relations (j2.39p . 

Let us rewrite the kernel '^>,"^p(2, ifja, /3) as follows (for brevity we will omit the normal- 
ization factor Am) 

{P^w-'za-')l^^I'^{z,a)IP{w,p) = (w/ p-' a-'z^Y^^I'^' {z,a)IP' {w, l3) . (B.5) 
The key point is that the factor Xm = ( (3^^ a^^Za ) depends only on the variables z~i^ , 

\ ^ J mm 

with j > m and on the variables Wkj with j < m. It can be shown as follows (we consider 
only z case): X^, depends on the variables {za)km with k = m + 1, . . . , N , which, due to the 
triangularity of the matrix Za, are given by linear combinations of the elements {z~^)kp with 
p > m. The matrix z~^ is determined by the equation 

z^^ a^^ = d~l^ z~^ , (B.6) 

which follows from Eq. (j2.44p . It is straightforward to see from the above equation that the 
matrix elements {z~^)kn depend on the elements z~^ with j > k only. Therefore, taking into 
account the commutation relations of Eq. ()2.39p one can transform (IB.4P to 

^^-'P' f^!rj'p'(^'^!^'^) fc'r'«"'^-)^ I'''{z,a)lP'{w,(3) 

= [c^'r'a-^aY^^n^'p' (7^^'"^,(r?,e|z,«))T-'(z,a)J'''(u;,/3)) 

where on the last step we use the transformation (jB.Sp . 
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C Appendix: Proof of the estimate (14.161) for matrix elements 
of factorizing operators 

Here we derive the estimate (|4.16p for the matrix element of the factorizing operator 'JZ^^{X). 
Let en{z) and ek{w) be the basis vectors ()3.5p in the quantum and auxihary spaces, respectively. 
Defining 

i>k 

one can represent en{z) as 

e„(z) = nkn{o-)Ek{da)I'^ {z, a)\o,=i(c,,j=o,i>j) , (C.2) 

where ^^^(cr) = i}fj(ek,en) and the sum over repeating (multi)indices is implied. The matrix 
element [n[l\x)]'!^;;'' can be expressed as follows 

[7^S?(ACL"' = nni{p)nrnj{a)Ern'{d,)En'{d^)Ej{d^)Ei{df^)ni^'hw,z\^^ 

' ^ z=w=a=(3=l 

(C.3) 

We are interested in the behavior of this matrix element at rijk — > oo, j = k + 1, . . . , N , all 
other indices mji,m'j^ and riji, n'j^, i ^ k being fixed. It is clear that the multi-index j in r.h.s. 

of ()C.3P varies in finite limits. Using Eq. (f3:22|l for 7^J^p and carrying out differentiation with 
respect to a and z one finds that the matrix element ()C.3p is given by a sum of terms which 
have the form 

nME^,{d^)Ei{d-p)CPw)ll'' P{w,B) , (C.4) 

w=f3=l 

where K is some constant. The polynomial P{w,P) has a finite degree and does not contain 
large factors rijk- The factor {I3w)kk = 1 + "^pPpkWpk depends only on the "large" variables 
ujpk, Ppk- After differentiation with respect to all other variables one gets for ()C.4p 

^MU^^—Aa;;-) CMlk"" Pin^pkJpk) ^ , (C.5) 

where again the polynomial P does not contains large factors njk- Finally, one gets that (jClK]) 
is given by the (finite) sum of the terms 

o f , r(-A + K + Y.f=k+Mk - Sj)) , / . 

[nk+ik - Sk+iV- ■ ■ ■ [n^i, - snV- 

where sj are some constants, i? is a polynomial and the difference of indices ijk — n'^^. is finite 
at n'-^^ oo, j = k + 1, . . . , N. 

Let us estimate the coefficient Qni{p) when all indices except rijk and ijk, j > k are small. 
To this end we will use the integral representation (|3.15p . 



^ni{p)=CN{p) / DznN{z)en{z)ei{z) , (C.7) 
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with /i AT (z) = n^i^ ^ 



Pj—pj+i. We recall that by definition the integral (IC.7P 
for arbitrary {Aj} is understood as an analytic continuation from the region of A's where it 
converges. To stress that the matrix element f2„j(p) is considered in the special "kinematic" 
{njk,jjk — ^ oo) we put more labels on it, Q.ni{p) — > and switch from p to A. 

To get the necessary estimate we proceed as follows: First we show that the coefficient 



J7^''^(A) for /c > 1 can be represented as the sum of integrals ()C.7p for sl{N — 1) case 



(C.8) 



It is important that the sum in (jC.SP over q goes in a finite range and the coefficients Cq and 
parameters {Xq)j = \j + Sj do not depend on njk,ijk- Continuing this procedure one can 
represent ^}^1^{X) as the sum of the elements O^'j^(Ag), where M = — A; + 1, thus reducing 

the problem to an evaluation of the matrix elements O^'-^(Ag). 

To prove (jC.SP let us make the change of variables z = w^^ and carry out the integrations 
over variables W21, . . . ,wni- Let us represent matrix w as 



w 



1 

a b 



( 1 

ai 

02 




1 

^21 





1 



0\ 







(C.9) 



VoAT-i 6iV-ll • • • bN^l,N-2 1/ 



and examine dependence of Ap on the elements ai, . . . ,aN-i. For N x N matrix M we define 
Ap(M) = det Mp and Mp is {N - p) x {N - p) matrix with elements Mij, i, j = p + 1, . . . , N. 
Noticing that Ap{z^ z) = Ap{ww^) = det \ {bb^)ij + aiaj\ij=p^,,,N-i one concludes that 



Ai = Ai(ai, . . . ,OAr_i) , A2 =A2(a2,...,aAr-i), ... AAr_i = Aat.i 
Let B = bb^ and Ap = (op, . . . , a7v-i)- Taking into account that 

Ap{ww'') = det(5p_i +Ap(g) Al) = Ap^i{bb^) ■ {1 + AI B-\ Ap 
one gets for the measure 



(oat- 



(C.IO) 



(C.ll) 



N-l 



PNiW) 



p=i 



N-l 



n A^Y {bb^){l + AlB;\Ap)-'r>-i 



p=i 



--PN-iib) Hil + AlB^l.Ap)-^"-' 
p=i 



(C.12) 



where fi^-iib) = Y\p=i ^p^''^^ ^i^b^) i'^oibb^) = !)• Now one can consequently carry out 



the integration over ai, 02, . . . , a^-i- Indeed, let us consider integral 

-Ai-l 



(C.13) 
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Representing i?.^.^ = niij/ detB, where rriij = {—iy~^^Mij, Mij being a minor of the matrix B, 
one derives 



1 + 4 ^1 ) = ( 1 + 
1 + 



1 



detS L 
mil 
detB 
mil 



I i2 t t t 

mil + a[mijaj + almnai + a-mijaj 



ai + mij/miiajf H \-al {mnmij - mumij} aj 

Tfl-i 1 



11 

,t^R-l> 



(C.14) 



Here the summation over repeated indices ( i,j = 2, . . . , N — 1 ) is imphed and we make use 
of the identity 



{mnmij — mumij} = mn det B ■ [B^ . 
After shifting and rescaling the integration variable one gets for ()C.13p 

min{m,rn) ^ ^ j_ n\ q+1 



(C.15) 



q=0 



det B Y 
mil J 



mijajy^-i f fnijUj 



detB 



detB 



m—q 



l + AlB^^A2)-^'+\ (C.16) 



where Cg{X) = C^C^q\T{X - q)/T{X + 1) and mu = det Si = Ai(66)t. Evidently, integrating 
over 02,03, •• • one again encounters the integrals of the type (IClSh . hence the final result of 
integrations can be cast into the form ()C.8p . This calculation also shows that the integral ()C.7p 
is a meromorphic function of Ai, . . . , Aat-i. 



In order to calculate the element ^l^-^{\q) one can repeatedly integrate over the variables 
^Nj, J > 1 in the last row, then over z^^ij, j > 1 in the row A'^ — 1 and so on. Representing 
the matrix z in the form 



/I ... 0\ 

621 1 ... 

^Af-ll bN-l2 ■ ■ ■ 1 

V oi 02 ... aN-i 1/ 



(C.17) 



{z^z)ij = {b^b)ij + OjOj, for i,j = 1, . . . , N — 1, one easily finds that Afc(z) depends only on 
the variables oi, . . . , o^. The integration over un-i, 0,^-2 • • • , 02 goes along the same lines as 
before. For instance, 



N-l „ N-1 

n / d-a, n A,T^-^(.t.) 

k=2 '' k=2 



IN -I 



'N-2 



k=2 



R{Z21,. . . ,ZNl) 



where 7Ar_i = vr^ ^(AAr_i(AAr_i + AAr_2) . . . (Aat^i + . . . + A2)) ^ and 

\t + |^;2l| + . . . + \ZN-ll\ ) 



R[Z21, 



.ZNl) 



[1 + |2;2lP + • • • + \ZN1 



|2\A2 + ...+Aj> 
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Therefore, one gets 



l<j<k<N-l 



N N-l / n \ / N \ -AiV-1 

— 1 r / -1 I \ ^ [ 1 2 1 111 \ ^ 1 ^ 1 2 



^pl ^pl 

q p=2 n=2 \ k=2 ) \ k=2 / 

where irip = ipi + rp{q), Ip = Upi + Sp{q), Hp = Xp + 6p{q) and Aat = /ii + . . . + /iAr_i. The 
parameters c{X,q),rp{q), Sp{q),5p{q) do not depends on nki,iki, k = 2,...,N and the over 
multi-index q goes in the finite hmits. FinaUy, after integration over Zpi, p = 2, . . . , N one gets 
for each term in the sum 

p=2 fc=3 r(i-Afc + Mfc) r(i-A2 + M2)' ^ ^ 

where A^. = /xi + . . . + and = rrik + . . . + niN- Together with Eq. ()C.6P this results in 
the estimate ()4.16p . 
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